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A new two-dimensional (2D) layered material, namely, titanium trisulfide (TiS 3 ) monolayer sheet, is pre¬ 
dicted to possess desired electronic properties for nanoelectronic applications. On basis of the first-principles 
calculations within the framework of density functional theory and deformation theory, we show that the TiS 3 
2D crystal is a direct gap semiconductor with a band gap of 1.06 eV and high carrier mobility. More remarkably, 
the in-plane electron mobility of the 2D TiS 3 is highly anisotropic, amounting to -10,000 cm 2 V -1 s _1 in the b 
direction, which is higher than that of the M 0 S 2 monolayer. Meanwhile, the hole mobility is about two orders 
of magnitude lower. We also find that bulk TiS 3 possesses lower cleavage energy than graphite, indicating high 
possibility of exfoliation for TiS 3 monolayers or multilayers. Both dynamical and thermal stability of the TiS 3 
monolayer is examined via phonon-spectrum calculation and Bom-Oppenheimer molecular dynamics simula¬ 
tion in NPT ensemble. The predicted novel electronic properties render the TiS 3 monolayer an attractive 2D 
material for applications in future nanoelectronics. 


The successful isolation of two-dimensional (2D) graphene 
in 2004 m has captivated great research interests in 2D mate¬ 
rials, particularly atomic-layered materials with weak inter¬ 
layer van der Waals bondingQ. Besides the grapheneQ]- 
0, the family of 2D layered materials also include transition 
metal dichalcogenides (TMDCs)|[2. |7J [8], hexagonal boron 
nitride (h-BN)[|9} jT0|, s i 1 i c ene BTT14131 . germanene lfl4ll . and 
phosphoreneil5j|T63, among others. These 2D atomic-layered 
materials exhibit not only novel 2D geometries as they repre¬ 
sent the thinnest crystalline solids that can be formed, but also 
new and exotic condensed matter phenomena that are absent 
in their bulk counterparts |2] [8). For example, single graphene 
sheet is a zero-gap semiconductor with a linear Dirac-like dis¬ 
persion near the Fermi level, while bulk graphite is known 
to exhibit semimetallic behavior with bandgap overlap of -41 
meV. 03 A M 0 S 2 monolayer sheet possesses a direct bandgap 
of - 1.8 eV, while bulk M 0 S 2 possesses an indirect bandgap of 
1.29 eVlfl8l. The bandgap of a few layer phosphorene (a 2D 
form of black phosphorus) is highly layer-dependent. For the 
phosphorene monolayer, the bandgap is ~1.5 eV while for the 
bulk phosphorus, the bandgap is merely -0. 3 evUmiEHZD. 

2D layered materials offer opportunities for a variety of ap¬ 
plications, particularly in next-generation electronic devices 
such as field-effect transistors (FET) and logic circuits. For 
high-performance FET applications, a 2D material should 
possess a moderate bandgap and reasonably high in-plane car¬ 
rier mobility. Graphene is a highly promising 2D material for 
high-speed nanotransistors due to its massless charge carriers. 
However, it lacks a bandgap for controllable operations. |[23l - 
[25l The molybdenum disulfide (M 0 S 2 ) monolayer sheets are 
more promising for FET applications since not only they pos¬ 
sesses a direct bandgap of ~ 1.8 eV. lfTSl but also the 2D M 0 S 2 - 
based FET devices show good performance with a high on/off 
ratio of -10 8 as well as a carrier mobility of -200 cm 2 V 's 1 . 
The latter can be enhanced even up to 500 crrrV 's 1 with 
improvement. 1130 Also, recent experiments demonstrated 
that FET devices built upon few-layer phosphorene exhibit 
reasonably high on/off ratio (up to 10 4 ) and appreciably high 


hole mobility of -55 cm 2 V 's 1 (at a thickness of -5 nm) 
to -1000 cm 2 V _1 s~ 1 (at a thickness of -10 nm).IT5l fl6) 
Nevertheless, new 2D layered materials with moderate direct 
bandgap and high carrier mobility are still highly sought. In 
this work, we show an ab initio calculation evidence of a new 
2D layered material - the TiS 3 monolayer sheet - with the 
desired electronic properties. 

Historically, bulk materials such as graphite, TMDCs and 
black phosphorous were studied well ahead of their 2D 
layered-material counterparts. Likewise, properties of bulk 
TiS 3 are known much earlier than those of 2D form. Bulk 
TiS 3 has a monoclinic crystalline structure (with the space 
group of p21/m), and the TiS 3 crystal can be viewed as stacked 
parallel sheets with each sheet being composed of ID chains 
of triangular TiS 3 unit. These sheets interact with one an¬ 
other via the van der Waals (vdW) forces. Il28l [291 It is also 
known that materials with stacking-layer structures can be 
a good precursor for contriving 2D atomic layers via either 
exfoliation P()l [31 1 or mechanical cleavage |32) . Several elec¬ 
trical and transport measurements have been reported, P3| [34) 
showing that the bulk TiS 3 is an n-type semiconductor with 
carrier mobility of -30 cm 2 V _1 s _1 at room temperature. The 
mobility can be further enhanced up to -100 cnrV^'s -1 at 
the low temperature 100 K. P4) Moreover, optical absorption 
measurements indicate that the bulk TiS 3 exhibits an opti¬ 
cal gap about 1 e V. 051 More importantly, several recent ex¬ 
periments demonstrate that thin films of TiS 3 with thickness 
of 10 2 nanometers possess a direct bandgap of -1.1 eV and 
exhibit good photo response. P6U38) The moderate bandgap 
of bulk TiS 3 coupled with relatively high carrier mobility 
renders the bulk TiS 3 a highly promising precursor for iso¬ 
lating 2D TiS 3 sheets with desired properties for nanoelec¬ 
tronic applications. For the 2D TiS 3 monolayer sheet, ge¬ 
ometrical optimization and electronic structure calculations 
are carried out using density-functional theory (DFT) meth¬ 
ods within the generalized gradient approximation (GGA) and 
with the Perdew-Burke-Emzerhof (PBE) exchange correla¬ 
tion functional, as implemented in the Vienna ab initio sim- 
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FIG. 1. (Color online)(a)A 2x2x1 supercell of the bulk TiSj struc¬ 
ture, (b) top view of a 2x2 T 1 S 3 monolayer sheet (left), and the first 
Brillouin zone and the high symmetry points associated with the 
monolayer (right). The grey and yellow spheres refer to Ti and S 
atoms, respectively. 


ulation package (VASP). ||39| The Grimme’s D2 dispersion 
correction [ ;401 is adopted to account for the long-range vdW 
interactions. The ion-electron interaction is treated using the 
projector-augment-wave (PAW) technique and a kinetic en¬ 
ergy cutoff of 500 eV is chosen. A vacuum space of ~20 
A along the direction normal to the monolayer plane is un¬ 
dertaken so that the interlayer interaction due to the periodic 
boundary condition can be neglected. For the geometric op¬ 
timization, a 7x10x1 Monkhorst-Pack RT) grid is used and 
all structures are relaxed until the forces on the atoms are less 
than 0.01 eV/A, and the total energy change becomes less than 
1.0 X 10 5 eV. For total energy calculations, a fine 35x50x1 
grid is adopted. Since the PBE functional tends to underes¬ 
timate the band gap of semiconductors, the hybrid HSE06 
functionalj42] is also used to compute the band gap of op¬ 
timized TiS 3 monolayer sheet. 

The carrier mobility (/ 1 ) is calculated based on the defor¬ 
mation theory proposed by Bardeen and Shockley. m Due to 
the fact that for inorganic semiconductors, the coherent wave¬ 
length of thermally activated electrons or holes is close to the 
acoustic phonon wavelength and is much longer than typi¬ 
cal bond length, the scattering of a thermal electron or hole 
is dominated by the electron-acoustic phonon coupling. |431 
The deformation theory has been widely used to evaluate p of 


low dimensional systems. il9ll44lj48l On the basis of effective 
mass approximation, the charge mobility in 2D materials can 
be expressed as: 


2 eh 3 C 
3k B T\m*\ 2 E\ 


(1) 


Here, C is the elastic modulus defined as C = [<9 2 £y<9(!> 2 ]/S 0 , 
where E is the total energy of the system (per supercell), and 
6 is the applied uniaxial strain, and S° is the area of the op¬ 
timized 2D structure, m * is the effective mass, which can be 
given as m* = h 2 (d 2 E/dk 2 )~ l (where h is the Plancks constant 
and k is magnitude of the wave-vector in momentum space), 
T is the temperature, and E\ is the deformation potential con¬ 
stant, which is proportional to the band edge shift induced by 
the strain. E\ is defined as 6E - E\ x ( 61/Iq ), where 6E is 
the energy shift of the band edge position with respect to the 
lattice dilation 61/1 0 along either direction a or b, the energies 
of the band edges are calculated with respect to the vacuum 
level. 

The PBE-D2 optimized structure of bulk TiS 3 is shown in 
FigJTJa), and the associated lattice constants are a = 4.982 A, 
b = 3.392 A and c = 8.887 A, and lattice angle£=97.24°, all in 
very good agreement with the experimental results, a = 4.958 
A, b = 3.401 Aand c = 8.778 A, and £=97.32°. |28) Further¬ 
more, the computed band structures of the bulk TiS 3 from both 
PBE-D2 and HSE06 are shown in Supplemental Material 11491 . 
PBE-D2 and HSE06 give qualitatively the same results except 
for the band gap. Both PBE-D2 and HSE06 calculations indi¬ 
cate that the bulk TiS 3 is an indirect gap semiconductor from 
r (0, 0, 0) to Z (0, 0, 0.5). The PBE-D2 computation gives a 
band gap of 0.21 eV, while HSE06 gives 1.02 eV. The latter 
agrees well with the measured optical gap which is around 1 
eV.(35) The good agreement between the benchmark calcula¬ 
tions and experiments for the bulk TiS 3 show that the theoreti¬ 
cal methods chosen for this system is reliable. In addition, the 
band structures near the conduction band minimum (CBM) or 
the valence band maximum (VBM) of the bulk TiS 3 exhibit 
notable in-plane dispersion behavior (from T to Y (0, 0.5, 0)) 
or F to B (0.5, 0, 0)), indicating that the 2D TiS 3 monolayer 
sheet may have relatively high carrier mobility. The computed 
HSE06 band structures and density of states (DOS) of the 
TiS3 monolayer are shown in Fig[2ja). Since the original Z 
point of the bulk TiS 3 folds back to the F point for the TiS 3 
monolayer, the TiS 3 undergoes an indirect-direct transforma¬ 
tion from an indirect band gap semiconductor for the bulk to 
a direct band gap semiconductor for the 2D monolayer coun¬ 
terpart, akin to the case of MoS->. |iT8l Both VBM and CBM 
are located at the F point, yielding a direct band gap of 1.06 
eV. Moreover, from the orbital and atom projected DOS, we 
can see that the valance bands exhibit strong hybridization be¬ 
tween the S p states and Ti d states from -2 eV to the top 
of valence band, while the conduction bands are mainly con¬ 
tributed from the d states of Ti (see Fig(2ja)). The isosurface 
plots of the VBM and CBM are shown in Fig|2] (b), which 
show that the holes (from VBM) favor the a direction, while 
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FIG. 2. (Color online)(a) Computed FISE06 band structures of TiS 3 
monolayer sheet; T (0.0, 0.0, 0.0), Y (0.0, 0.5, 0.0), A (0.5, -0.5, 0.0), 
B (0.5, 0.0, 0.0) refer to the special points in the first Brillouin zone; 
red circles, green rectangles refer to the contributions of Ti d states 
and S p states; and the Fermi level is set to zero; (b) iso-surface plots 
of the charge density of VBM (left) and CBM (right) of the TiS 3 
monolayer sheet, with an iso-value of 0.003 e/Bohr 3 . 


the electrons (from CBM) favor the b direction. To compute 
the 2D elastic modulus (C) and the deformation-potential con¬ 
stant (Ei), we dilate the lattice of the cell up to 1.5% along 
both a and b directions, and then calculate the total energy 
and the positions of CBM and VBM with respect to the dila¬ 
tion. The atomic positions are relaxed at the dilation, and the 
electronic energies are calculated at the PBE-D2 level with 
the ultra-fine k-meshes (35x50x1). We note that although 
the PBE functional underestimates the band gap, it can give 
quite good carrier mobility data for M 0 S 2 , m graphene, 1461 
graphyne{45j and graphdivne l47l . The total energy-strain re¬ 
lation and the positions of CBM and VBM with respect to 
the strain are plotted in Fig|3] As shown in Fig|3]b), the re¬ 
sponse of CBM and VBM to the applied strain appears to be 
highly anisotropic. The CBM increases monotonously with 
the strain either along a or b direction, while the VBM de¬ 
creases monotonously with the strain along the b direction but 
increases along a, resulting in band gap increase due to the 
strain along b but band gap decrease along a. The 2D modu¬ 
lus (C) is attained by the quadratic fitting of the total energy 
versus strain, and the deformation potential constant (E\) is 
calculated by the linear fitting of the CBM (VBM)-strain re¬ 
lation. With C, £1 and the effective mass known, the carrier 
mobilities are calculated via equation (1). These data and the 
relaxation time (r = pm*/e) are summarized in Table I. 
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FIG. 3. (Color online)(a) Strain-total energy relations and (b) shifts 
of VBM and CBM under uniaxial strain along a and b directions for 
TiS 3 monolayer sheet, 61 refers to the dilation along a or b, while Iq 
refers to the lattice constant of a or b at equilibrium geometry. In (b), 
the vacuum level is set at zero for reference. 


TABLE I. Calculated deformation-potential constant (£1), 2D mod¬ 
ulus (C), effective mass (m*). relaxation time (r), and electron and 
hole mobility (p ) in a and b directions of TiS 3 monolayer sheet at 
300 K. 



Ei (eV) 

C (N/m) 

m * (m e ) 

r (ps) 

p (10 3 cnrWV 1 ) 

electron (a) 

0.73 

81.29 

1.47 

0.84 

1.01 

hole (a) 

3.05 

81.29 

0.32 

0.22 

1.21 

electron ( b ) 

0.94 

145.05 

0.41 

3.23 

13.87 

hole ( b ) 

-3.76 

145.05 

0.98 

0.085 

0.15 


As shown in Tahiti] the 2D modulus along b is nearly two 
times higher than that along a direction. This is because the 
Ti-S bond strength along a is weaker than that along b direc¬ 
tion as the Ti-S bond length is 2.65 A along a while 2.45 A 
along b. The difference between the Ti-S bond strength along 
a and b also makes the deformation-potential constant along b 
larger than that along a, as the band energies are more sensi¬ 
tive to dilations along b than along a. The effective mass also 
shows an anisotropic feature: in a direction, the effective mass 
of hole is much smaller than that of electron, while in b direc¬ 
tion, the effective mass of electron is about half of that of the 
hole. These results can be well explained by a charge-density 
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plot of CBM and VBM as shown in Figure 2(b) where one can 
see that VBM electrons are quite localized along b while the 
CBM ones are delocalized along b but localized along a. 

The predicted carrier mobilities of TiS 3 monolayer are 
highly anisotropic. The computed electron mobility along 
b direction is 13.87xl0 3 crrrV 's ', about 14 times higher 
than that along a direction (l.OlxlO 3 cnrV^'s -1 ), while 
the hole mobility along a direction is 1.21xl0 3 cnrV 's~\ 
about 8 times higher than that along b direction (0.15xl0 3 
cnrV^'s -1 ). It is worthy of mentioning that the pre¬ 
dicted carrier mobilities are notably higher than those of the 
M 0 S 2 monolayer sheet (which are in the range of 60-200 
crrrV 's 1 ) ||44| |. Especially, along the b direction, the elec¬ 
tron mobility is about 100 times higher than the hole mobility, 
making b direction more favorable for the electron conduc¬ 
tion. The large difference in electron/hole mobility can be 
exploited for electron/hole separation. 

Although TiS 3 monolayer exhibits some novel properties 
for potential nanoelectronic applications, feasibility of isola¬ 
tion of TiS 3 monolayer sheet via either exfoliation or mechan¬ 
ical cleavage technique has yet to be confirmed. To example 
this feasibility, we calculate the cleavage energy by introduc¬ 
ing a fracture in the bulk TiS 3 (see Fig |4jb)). To this end, 
the total energies under variation of the separation d between 
the fractured parts are computed to simulate the exfoliation 
process. 150 1511 The resulting cleavage energy is plotted in 
Fig|4](a). It can be seen that the total energy increases with 
the separation d and gradually converges to the ideal cleav¬ 
age cohesion energy of about 0.20 J/m 2 . The latter value is 
less than the experimentally estimated cleavage energy for the 
graphite ( ~0.36 J/m 2 ) j52l , indicating that the exfoliation of 
bulk TiS 3 should be highly feasible experimentally. The sta¬ 
bility of TiS 3 monolayer is another issue that should be ex¬ 
amined. First, we compute the phonon spectrum of the TiS 3 
monolayer, based on density functional perturbation theory 
with the linear response as implemented in the QUANTUM- 
ESPRESSO package.J53l As shown in Figj4] (c), the TiS 3 
shows no imaginary phonon mode, indicating its dynamical 
stability. Next, we perform BOMD simulations. The constant- 
temperature (300 K) and -pressure (0 GPa) ( NPT ) ensemble 
is adopted. Here, the time step is 2 fs and the total simulation 
time is 8 ps. As shown in Fig|4](d), the in-plane structure in¬ 
tegrity of TiS 3 monolayer is well kept during the BOMD run, 
suggesting good thermal stability of the TiS 3 monolayer. 

In conclusion, we predict a new 2D material, TiS 3 mono- 
layer sheet, which is a semiconductor with a desired direct 
band gap of ~1 eV. The electron mobilities of the TiS 3 mono- 
layer are dominant and highly anisotropic. More specifically, 
the electron mobility along b direction exhibits a very high 
value of 13.87xl0 3 crrrV 's ', rendering the T’iS 3 monolayer 
particularly attractive for future applications in nanoelectron¬ 
ics as its mobility is even notably higher than that of the M 0 S 2 
monolayer. The computed ideal cleavage cohesion energy 
for TiS 3 is about 0.20 J/m 2 , less than that of graphite (-0.36 
J/m 2 ), indicating the isolation of a 2D TiS 3 monolayer can be 
technically attainable via either liquid exfoliation or mechan- 
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FIG. 4. (Color online)(a) Cleavage energy E c; as a function of the 
separation between two TiS 3 monolayers, (b) Schematic view of the 
exfoliation of a TiS 3 monolayer from bulk, (c) Phonon band structure 
and density of states of TiS 3 monolayer and (d) Snapshot of TiS 3 
monolayer at 8 ps of the BOMD simulation in the NPT ensemble at 
300 K and 0 GPa. 


ical cleavage as done for isolation of 2D graphene or M 0 S 2 
sheet. 13014321 Lastly, dynamic and thermal stability of TiS 3 is 
confirmed by both phonon spectrum and BOMD simulations. 
Thus, we expect that fabrication of 2D TiS 3 monolayer and 
measurement of its electronic properties will be likely accom¬ 
plished in the near future. 
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